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Abstract 

The consistency of effective models with QCD is investigated through the 

use of the QCD sum rule. Taking the potential model for the heavy quark 

system, we apply the method to two phenomenologically successful parameter 

sets, and obtain the dependences of the model parameters on the QCD scale 

A. Comparison with the expected scaling laws allows us to reject one of the 

two sets. The method is applicable to any model which reproduces the low 

lying spectra of hadronic systems. 
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Heavy quark systems have been the most successful place in the study of strong interac- 
tion physics. There exist many methods that give quantitative agreement with experiment, 
both in the QCD based methods and in the QCD inspired models. As for the former, the 
heavy quark effective theory has related various observables in the framework of perturbative 
QCD, and the QCD sum rule has been applied to various channels to calculate the masses of 
the heavy mesons. As for the latter, the non-relativistic quark model(QM) has proved itself 
to be accurate enough to reproduce a vast number of physical observables, with the bonus 
that it stays very close to intuition. It is the mass of the heavy quark itself that sets the 
scale in the meson masses, the role of the strong interaction being only to give the splitting 
among them. 

A problem about the QM is that there are usually several sets that seem to work equally 
well in a phenomenological sense. Being a model, it is natural to expect that it reflects 
the effective degrees of freedom when all the complicated degrees of freedom in QCD are 
integrated out, and ultimately its parameters be calculated starting from the fundamental 
theory. At this stage, no such calculation has been done, which leaves us to rely either 
on their self-consistencies [^, or on their phenomenological successes. In this paper, we 
propose the use of QCD sum rules to investigate the consistency of various models with 
QCD. Evaluation of the vector current polarization function in the cc channel enables us 
to exclude one of the two specific sets in the potential model for the charmonium system. 
The method is in principle applicable for any phenomenological model, the only requirement 
being that it reproduces the observed mass spectrum and leptonic decay widths. 

The QCD sum rule basically deals with the polarization function, evaluating it on the 
one hand by the operator product expansion(OPE) in perturbative QCD and by a phe- 
nomenological model representing the observed spectrum on the other. Here we consider 
the vector polarization function 

/ + 00 
d'x e'^-^ (0| T {j, (x) J, (0)) |0) = {q,q, - g V)n(g') (1) 

where j^(a;) = '?/'c(a;)7^'?/'c(a;) and is the 4-momentum carried by the system. Instead 
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of extracting the physical masses by using the simphfied spectral function, we substitute a 
more realistic spectrum calculated by the model on the phenomenological side to see how the 
parameters in the model vary with the quantities appearing in the OPE. More specifically, 
we calculate the variation of Il{Q^) with the QCD scale parameter A and fit it for a certain 
range of with the variations in the QM side with respect to the parameters of the model. 
It is essential that we fit the change in the polarization function and not the function itself, 
since the latter is usually saturated by the first few resonances, and is insensitive to the 
parameters as long as they accurately reproduce the first few spectra. The validity of the 
model parameters is investigated through their A dependences. 
For our analysis, we take the simple 'Cornell' type potential 

V{r) = -- + Kr + Vo, (2) 
r 

with two sets of parameters as our candidates. Set A has four adjustable parameters and 
reproduces m^(i5),m^(25), r^(i5)^e+e- and r^(2S')^e+e- , with the radiative correction factor 
(1 — ^a) for the widths The parameters take the values 

{m = 2.04[GeV], a = 0.579[1], k = 0.172[GeV^], Vq = -1.12[GeV]}, 

which are essentially the ones employed by Set B, on the other hand, has one restriction 
a = = 0.27, and is thus adjusted to reproduce only mis,m2s and r^(i5)^e+e-- The 
parameters are given by 

{m, a, K, Vo} = {1.78, 0.270, 0.222, -1.00} 

. Their characteristics are summarised in Table |. 

The former is naturally better in reproducing the experimental data, and is the most 
frequently used set in modern quark model analyses. The latter on the other hand has the 
strong point that the Coulomb force is directly correlated with the one-gluon exchange of 



*A stands for Amom in this particular case. 



2 



QCD. Set B with the couphng with DD channel will be discussed later. With these models 
at hand, we calculate 

^ ' 

with Imn(s) given by Imn(s) = E f K0| j^(0) |0,n)|^5(s - M^). \k,n) denotes the 
nth bound state with momentum k, with the normalization given by (p, n \ p', n') = 
{27i)^2pq6^ {p — p')6nn'- ^ iucludes the possible polarizations of the state. The matrix el- 



ements are evaluated with the use of the van Royen-Weisskopf formula [|TT 



On the other hand, the OPE side of the QCD sum rule is calculated in the usual manner 

V 

n(g^) = CoiQ'') + CaiQ'') (0| -G,,G,, |0) + .... 

where the dots indicate higher dimensional contributions Neglecting light quarks, we have 
two parameters on this side, i.e. the QCD scale parameter A and the heavy quark current 
mass, rric- We are here interested in the variants of 11 with respect to A, with fixed. 
The essential point is that we control the strength of the interaction solely through the 
parameter A. This parameter comes in through the quantities as and the gluon condensates 
( (0| ^G^yG^j_y |0) etc.). (0| ^Gfj_yG^y |0) behaves as A^ since it has no anomalous dimension 
(we are neglecting light quarks) and as will behave as l/ln(mc/A). It is easy to see that 
taking A ^ will give for both quantities, which means that the interaction is turned off 
and we are left with only the free particle contribution to the polarization function. 
Now we can compare both sides concerning their variants with respect to A, 

gn(g^) _ gn(g^) dm gn(g^) 5a 5n(g=^) , 5^ 
~^r~ ''"^ ~ '''•'^ 9A '"^^ + dk '"^ + l'^''" dk ''"^ ■ 

The derivatives of n(g^) with respect to the parameters are calculated by numerical differ- 
entiation, at some fixed values of A, mc,a, etc. The coefficients ( ^ \m.ci--- ) ^ire obtained 



''^We take into account upto operators of dimension 8. Heavy quark condensates are reduced to 
gluon condensates with the use of the heavy quark mass expansion ||2|,Q . 
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as a result of the fit in the region = 0.5 ~ 3.0GeV^ 0. In the actual calculation, we 
evaluate the moments of n(Q^), i.e. n„(Q^) = ^{—-^)^Il{Q'^) since this has an effect 
of emphasizing the low energy part of the physical spectrum [^. We also evaluate for the 
parameter O, 

instead of for reasons of convenience. This implies that the parameter behaves as O ~ A'*" 
near the actual value of A. 

The fit is obtained in the simplest way through a linear fit formula (matrix inversion). 
We have found that the matrix is practically singular, which means that the effect of the 
parameters are redundant when reflected in the behavior of the polarization function. This is 
not surprising when we recall that the function itself is saturated by the first few resonances, 
which leaves us with only 2n degrees of freedom where n is the number of resonances required 
to saturate the function. Actually, the lowest one is dominant and we found that only two 
out of four parameters were independent. This forces us to fix two of the A dependences 
of the parameters, which we choose to be k ~ A^ and Vq ~ A. This corresponds to taking 
the quenched approximation for the interaction, which is in agreement with the spirit of the 
quark model |. 

For set A, fitting the curves shown in figure gives the scaling behaviors of the model 
parameters expressed as O ~ A'^o, with do 

{d^, da, 4, d^} = {0.28 ± 0.02, 0.6 ± 0.3, 2., 1.} 

where d^, and d^ are inputs. The error is due to the following ambiguities 

• Region of where the fit is performed. 



■'■Of course we do not explicitly handle the dynamical (radiative) effects of the gluons |13|. We 
take the view that their effects are integrated out and arise as a change in the potential and the 
quark mass when we restrict ourselves to the two fermion sector. 



• Which moment we take. : Hq ~ IIio 

• Value of ttsi'lml), or equivalently, Amom : 130 ~ 220 [Me V] 

• Value of (0| ^G^^G^^ |0) : (360 ± 50[MeV])^ 

Although we know of no correct formula for the analytic behavior of m against A, we can 
give an estimate for the A dependence from the following arguments. First, it is natural to 
expect that m reduces to rric (the 'on shell' mass ^ ) in the limit A ^ 0, which is the free 
limit. This determines the sign of dm- Second, from dimensonal analysis we expect that 
Am = m — rric ~ c^A where Cm is some constant of proportionality. This requires that 

d(lnm) Am 
dm = ry = 6 

which, after substituiton of m = 2.[GeV] and rric = 1.5 [GeV], gives dm = 0.25 which is 
clearly close to the value of set A. (allowing a dependence of m ~ ln(A/mc) does not alter 
the essence of this estimate.) Also the dependence of a on A is consistent with that of the 
QCD as parameter ((ia ~ daj- 

On the other hand, for set B we obtain 

{dm, da, 4, 4} = {-0.2 ± 0.1, -4.9 ± 1.9, 2., 1.} 

We notice that they fulfil neither of the requirements discussed above. They are not appro- 
priate firstly since their signs are the opposite, and secondly because we would not expect 
such a strong dependence of a on A. 

These results clearly show that we should take set A as our phenomenological model in 
order to meet the requirements of QCD. This is rather contrary to our intuition since it 
seems that set B is more easily justified from the viewpoint of perturbative QCD. We have 
considered the following facts that might be the reasons for the failure of set B: 

1. Set B does not reproduce r^(25)^e+e- (Table |). The effect of the second bound state 
accounts for a non-negligible portion of the polarization function. 
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2. The 'radiative' correction (1 — ^a) differs largely in the two cases. Taking a to be the 
value of set A, this gives a correction of more than 50% of the magnitude. While this 
is often considered as a weak point in a quark model analysis, we on the other hand 
have obtained a result that supports this factor. 

To investigate the first effect, we have performed a coupled channel analysis similar to the 
one in For simplicity, we have taken only the DD channel into account, with the coupling 
potential 

{cc\V\DD) = g e-^^'^e,, 

where denotes the polarization vector of the cc system. This allows us to reproduce the 
decay width r^(25)^e+e-- Going through the same analysis but now with a larger parame- 
ter space {m, a, k, Vq, /i, g, we obtain figure ^ which shows the contribution from the 
changes in the model parameters fi,g,mD and k respectively. The difference in their mag- 
nitudes is clear. Requiring fi,g and niD to take appropriate values (ci^ = dg = l,dmo — -2) 
I , they clearly have only a small effect on the fit. We gain the values 

{d^, da, 4, 4} = {-0.15 ± 0.06, -3.8 ± 0.7, 2., 1.}, 

which is essentially unchanged from the previous fit. Thus, the modified set B (including the 
coupled channel effect) works equally well phenomenologically as the set A but the scaling 
behavior of the parameters is still inconsistent with that expected from QCD. 

One can easily diminish the effect of the factor (1 — ^a) by simply replacing Ilg with the 
ratio Un+i/Hn- The previous result obviously satisfies the fit, but due to the cancellation 
there might still be a set of parameters that is consistent with OPE. Notice that we have even 
fewer degrees of freedom (namely one), since taking a ratio means losing the information on 
the magnitude. The analysis gives the result 

{d^, da, 4, d,} = {0.10 ± .04, 0.2, 2., 1.}, 



dmjj is evaluated as {rriD — {rric + mu)}/m£i = 0.2. 
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where we have substituted da = 0.2. Varying da gives a change in the result which is included 
in the error. We have thus obtained the scaling behavior of set B close to that of set A, but 
this is only due to the fact that we are looking at a fit with fewer requirements. Requiring the 
individual moments to behave correctly will force the parameters to go to the inappropriate 
value. This indicates that it is indeed the 'radiative correction' factor (1 — ^a) that prohibits 
the use of set B. In other words, the factor is essential in achieving the consistency with 
OPE. 

In this paper, we have shown that it is set A of the model parameters that is consistent 
with the OPE, and is therefore appropriate for use in calculating physical quantities. This 
is in agreement with the choice of who justifies this set through comparison with exper- 
imental data. We have also shown that it is the 'radiative correction' factor itself that was 
responsible for this result, although its 'physical' meaning remains unclear. Our analysis 
shows that in one way or the other, the total magnitude of Imn(s) must be appreciably 
modified from the simple potential picture. 

The method is applicable in exactly the same fashion for a model that reproduces the 
spectrum and the leptonic decay widths. Applying finite temperature QCD sum rules, the 
method can also be used to extract the temperature dependences of the model parameters. 
In this case, one has to substitute two of the temperature dependences. 

T.M. wishes to thank N. Ishii for various ideas and discussions, S. Takeuchi for continu- 
ous encouragement, and finally K.Yamashita and T.Mizusaki for tips on numerical computa- 
tion. This work was supported in part by Grant-in-Aid for General Scientific Research (No. 
04804012) and by Grant-in-Aid for Scientific Research on Priority Areas (No. 05243102) 
from the Ministry of Education, Science and Culture. 
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TABLES 
TABLE I. Calculated Spectrum 





set A 


set B 


B with DD 


Exp. 


mis [GeV] 


3.097 


3.097 


3.097 


3.097 


m2S [GeV] 


3.685 


3.685 


3.685 


3.685 


miD [GeV] 


3.82 


3.775 


3.75 


3.770 


m404o[GeV] 


4.09 


4.14 


4.2 


4.040 


rv,(i5)-^e+e- [keV] 


5.26 


5.26 


5.26 


5.26 


^jP{2S)^e+e- [keV] 


2.14 


3.0 


2.14 


2.14 
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FIGURES 



Parameter Dependence of n(Q ) 
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FIG. 1. 

Variation with respect to A,m,a,K and Vq. Plotted are dlnUopE/dlnA, dlnHqM/dlnm and its 

counterparts. 
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FIG. 2. Variation with respect to K,fi,g and mo- Plotted are dlBllQM/dlnK,dlnIlQM/dlnfi 
and its counterparts. 



